Schiff moment of the Mercury nucleus and the proton dipole moment 
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We calculated the contribution of internal nucleon electric dipole moments to the Schiff moment 
of 199 Hg. The contribution of the proton electric dipole moment was obtained via core polarization 
effects that were treated in the framework of random phase approximation with effective residual 
forces. We derived a new upper bound \d p \ < 5.4 x 10 -24 e-cm of the proton electric dipole moment. 
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I. INTRODUCTION 

The interest to electric dipole moments (EDM) of ele- 
mentary particles and more complex systems like nuclei 
and atoms exists since 1950, when it was first suggested 
that there was no experimental evidence for symmetry of 
nuclear forces under parity transformation 1] . The inter- 
est was renewed after 1964 when it was discovered that 
the invariance under CP transformation, which combines 
charge conjugation with parity, is violated in if-meson 
decays. This provided a new incentive for EDM searches. 
Since the combined CPT transformation is expected to 
leave a system invariant, breakdown of CP invariance 
should be accompanied by T violation. Thus, there is a 
reason to expect that P and T violating EDMs should 
exist at some level. 

The experimental upper limit on the neutron EDM is 

| 

d n < 0.63 x 1(T 25 ecm. (1) 
The measured value for the proton EDM is Q 



d p = (-4 ±6) x 10~ M ecm, 



(2) 



and is compatible with zero. This corresponds to an up- 
per limit which is three order of magnitude weaker than 
the one for the neutron. 

The best upper limit on EDM ever obtained was in 
an atomic experiment with 199 Hg 3j. The result for the 
dipole moment of this atom is 



d( 19y Hg) < 2.1 x 1Q- 2 * ecm 



(3) 



Unfortunately, the implications of this result are some- 
what less impressive, due to the electrostatic screening 
of the nuclear EDM in this essentially Coulomb system. 
The point is that in a stationary state of such a sys- 
tem, the total electric field acting on each particle must 
vanish. Thus, an internal rearrangement of the system's 
constituents gives rise to an internal field E int that ex- 
actly cancels ~E ext at each charged particle; the external 
field is effectively switched off, and an EDM feels nothing 

Still, some P and T odd component of the electrostatic 
potential survives due to final nuclear size. It is created 



by the next moment in the nuclear electric dipole density 
distribution. This is the Schiff moment denned as Ifj 



(4) 



The Schiff moment generates a P and T odd electrostatic 
potential in the form 



(r) = 4ttS • Vtf(r) 



(5) 



Interaction of atomic electrons with the potential given 
by Eq. JSJ produces an atomic dipole moment 
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(6) 

Because of the contact origin of the potential, only the 
electrons in s and p atomic orbitals contribute to the 
dipole moment given by Eq. |JB}. 

The Eq. J3J is valid for any system of point-like charges 
e q . Let us split the sum in Eq.Q) into the sum over 
coordinates of nuclcons and the sum over coordinates of 
charges inside the nuclcons: 
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Here rjv is a nucleon position and p i is the position of 
the ith charge inside the nucleon. Combining the terms 
of the zeroth and first order in p and using J^. e^ = ejv, 
~Ylii e iPi = we obtain an expression for the Schiff 
moment as a sum of two terms. The first of them is 
similar to Q 
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(8) 



where ejv is equal to e for a proton and zero for a neu- 
tron. The mean value of this operator is nonzero only in 
the presence of the parity- and time-invariance violating 
nucleon-nucleon interaction 

The second term in the operator J7J) is related to the 
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internal dipole moments of the nucleons 



B. Core polarization 



S 2 = i]Td N (r%-(r% h ) 
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(9) 



The effects of the core polarization for a single particle 
operator can be treated by introducing a renormalized 
operator S satisfying the equation 



(i/u'lFfrv), (11) 



The previous calculations of the Schiff moment of a 
heavy nucleus 0, Q were performed in a simplified man- 
ner, without taking into account the residual interaction 
between a valence nucleon and the core nucleons. Only 
recently more microscopic studies of the Schiff moment of 
199 jjg Q an( j 225 p^ a a pp eare( j where effects of the core 
polarization with the effective forces for 199 Hg and the 
octupole deformation for 225 Ra based on Skyrm-Hartree- 
Fock method were discussed. In this work we would like 
to concentrate on the nucleon EDM contribution to the 
Schiff moment of the 199 Hg nucleus. In the picture of in- 
dependent particle model only an EDM of a valence nu- 
cleon contributes to the Schiff moment. In case of 199 Hg 
it is a neutron EDM. However, when a residual quasi- 
particle interaction between the valence neutron and the 
protons in the core is taken into consideration, the proton 
EDM contribution to the nuclear Schiff moment becomes 
nonzero. We calculated this contribution using a random 
phase approximation with effective forces. From the re- 
lation between the Schiff moment and the electric dipole 
moment of the Hg atom [Toj the new upper limit on the 
proton EDM was obtained. 



II. OUTLINE OF THE THEORY 

A. Nuclear mean field 

In our calculations we used full single-particle spec- 
trum including continuum. The single-particle basis was 
obtained using partially self-consistent mean-field poten- 
tial of The potential includes four terms. The 
isoscalar term is the standard Woods-Saxon potential 



Mr) 



V 



1 + cxp^ 



(10) 



with the parameters V = 52.03 MeV, R = 1.2709A 1 / 3 
fm, and a = 0.742 fm. Two other terms Ui s (r), and U T (r) 
were obtained in a self-consistent way using two-body 
Landau-Migdal-type interaction of [12$ for the spin-orbit 
and isovector parts of the potential. The last term is the 
Coulomb potential of a uniformly charged sphere with 
R c = 1.18A 1 / 3 fm. The mean field potential obtained in 
this way produces a good fit for single particle energies 
and r.m.s. radii for nuclei in the region around 208 Pb. 



where S is the bare Schiff moment operator given by 
JHl or • rift and e M are single particle occupation num- 
bers and energies. For static moments the external fre- 
quency uj — > 0. The value of the Schiff moment is given 
by the diagonal matrix element of the z-component of the 
renormalized operator (|11|) between mean field states of 
the last unpaired nucleon with a maximal angular mo- 
mentum projection. 



S = (/ij m = j\S z \[xj m = j). 



(12) 



For the residual interaction F we use the phenomeno- 
logical Landau-Migdal interaction that has the form 

F = C [g s {<Tx ■ <x 2 ) + g' s {a x ■ £r 2 )(ri • r 2 )] S(r 1 - r 2 ), 

(13) 

where C = 300 MeV fm 3 . The values of the empirical 
interaction constants g s and g' s are crucial for our calcu- 
lations. The proton contribution is proportional to the 
proton - neutron interaction g s — g' s . The constant g' s is 
determined from magnetic properties of nuclei and po- 
sitions of Gamow- Teller resonances. Its adopted value 
varies between g' s = 0.9-1 de pending on details of the 
mean field potential used [Tj, Uj, HI] • The constant g s 
is not so well defined. The magnetic moments and Ml 
transitions are to large extent isovector and they do not 
fix g s . An attempt to fix it from the structure of high 
spin states in 208 Pb has been done in ^(|. They found 
that g s = 0.25 had to be used in order to reproduce the 
excitation energies of the M12 and MIA states. Another 
value g s — 0.19 was quoted in the review paper |l5j . 

The Schiff moment operator can be presented in coor- 
dinate space in the form 



(i) 

m '■ 



(14) 



i=l 



where we have introduced the set of linear independent 
tensor operators 



T 
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Y^/(n), T 



(2) 
JM 



Y^(n), (15) 



where Yj M (n) is the vector spherical harmonic. 



For 



J = 1 we have T-, 



(i) 



-.(2) 



and T lm 



i(n-o-) 



-\ m '~ cr m , and ~ n m (n ■ erj - |cr TO . 
For a spherical nucleus we can separate the angular vari- 
ables and solve the obtained equations in coordinate 
space. The equations are 



S ai ( r ) = SS\r) + / A aibj (r,r')S bj (r')dr', (16) 
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Fig. 1: Core polarization effects in the neutron Schiff moment 
operator. The solid curve is S" ll (r); the dashed curve is the 
bare operator S'o 1 (r). 



Table I: Values of s p and s n for different g s and g' g . 
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dependence of the proton contributions induced by the 
core polarization. Note 1 order of magnitude difference 
in the scales in Fig. 1 and 2. The full curve in Fig. 2 
is the radial dependence at the first operator and 
the dashed curve is the radial dependence at the second 
operator Tj®. S pl (r) changes sign inside the nucleus, 
therefore its mean value is smaller than the mean value 
of S p2 (r) which is mostly negative inside the nucleus. 
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Fig. 2: Radial dependence of the proton effective Schiff mo- 
ment operators induced by the core polarization. The solid 
curve is S pl {r); the dashed curve is S p2 (r). 

where a = p,n and S§ % (r) is the radial part of the Schiff 
moment operator 10 multiplied by r. The kernel of the 
integral equation A ai i{r, r') was calculated by means of 
the Green functions of the radial Schrodinger equation. 



xr'R b K (r') (G>,r'|e K + u) + G*,(r, r'\e K - u )) , (17) 

where g™ = g nn = g s + g' s , g pn = g n P = g s - g' s , R b K {r) 
are the radial wave functions, and n b K are the occupation 
numbers. 

The solutions of Eq. (fTT|> for S al {r) are shown in Fig.l 
and 2. Figure 1 demonstrates the magnitude of the 
core polarization effects. Repulsive residual interaction 
(|13|l leads to not very significant decrease of the mean 
value of the Schiff moment. Figure 2 shows the radial 



III. RESULTS 

The value of the Schiff moment of 199 Hg can be pre- 
sented as a sum of proton and neutron contributions 

S — Spdp -\- Sjidji- 

In Table I we list the values s p and s n calculated for dif- 
ferent combinations of g s and g' s . From Table [I] one can 
see that the uncertainties in s p and s n due to uncertain- 
ties in g s and g' s are 

s p = 0.20 ±0.02 fm 2 , s n = 1.895 ± 0.035 fm 2 . (18) 

The main contribution to s p and s n comes from the sec- 
ond term in Eq. ©. The contribution of the first term 
is only -0.7 fm 2 in s n and 0.006 fm 2 in s p . 

The constraint for the Schiff moment of Mercury nu- 
cleus from the experiment Q can be obtained using the 
results of Ref. 10]. They calculated EDM of an atom 
created by the nuclear Schiff moment. For 199 Hg they 
found 

d= -2.8 x lO" 17 ^ [efm 3 ]. 

From Eq. we obtain the following upper bound for the 
Schiff moment 

IS^Hg)! < 0.75 x 10~ n e fm 3 , (19) 

From Eq. 1|19H we can give the following constraints for 
EDM of nuclcons: 

|dp| < 3.8 x 10~ 24 e cm, |d„| < 4.0 x 10" 25 e cm. (20) 

The constraint for the neutron EDM is worse than the 
existing result d n < 0.63 x 10~ 25 e cm 0, therefore, we 
shall not discuss it below. For proton EDM the esti- 
mate l|20f) is 1 order of magnitude lower than the exist- 
ing experiment d p = (—4 ± 6) x 10~ 23 e cm 0- In these 
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circumstances the question about a real theoretical ac- 
curacy of our approach becomes important. It is clear 
that the value ±0.02 cited in Eg. I|18fl does not reflects 
the real accuracy of the theory. It just came from the 
difference in adopted values of g s and g' s . The theoret- 
ical uncertainty appears from two sources. First, it is 
an uncertainty in the atomic calculations that couple the 
nuclear Schiff moment and EDM of an atom. We shall 
not discuss it here referring to the work [ljj. Second, it 
is an uncertainty in calculations of the core polarization 
effects using RPA with effective forces. The latter can 
be estimated from the following considerations. Using 
RPA with the effective forces we can fit different nuclear 
moments in one nucleus. Then, in neighbor nuclei the 
calculated moments will differ from the data. This dif- 
ference can be regarded as an uncertainty in the theory. 
In our experience this difference is of the order of 20% on 
the average, reaching sometimes the value of 30% [l7j . 
To be safe, we can adopt a conservative 30% uncertainty 
in calculations of s p . Therefore, instead of (|18f) we would 
prefer to write for s p 

s p = 0.2 ±0.06 fm 2 . (21) 



Since the error in i|21|) is not statistical, we can not give 
a probability distribution for s p . If one takes 0.14 fm 2 as 
a minimal value of s p , then it gives the following value 
for the proton EDM upper bound 

\d p \ < 5.4 x 10~ 24 e-cm. (22) 

In summary, we calculated the contributions of the pro- 
ton and neutron EDM to the Schiff moment of 199 Hg. 
The effects of core polarization were accounted for in the 
scope of RPA with the effective residual forces. A new 
upper bound of the proton EDM has been obtained from 
the upper bound on the atomic EDM of 199 Hg atom. 
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